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PLN RO DUCT LON 


Whether one 1s extracting information, such as separating 
Signals from noise in a communications system, detecting 
targets in a radar system, or preparting signals for trans- 
mission or display, Signal processing is a major portion 
of system design. This processing may be accomplished by 
several means, ranging from simple RLC circuits to complex 
digital computers. Several factors will influence the 
choice of processing scheme, among these are nature of the 
Signal and noise, desired results, weight, power consumption, 
cost, and convenience. Each scheme has its own advantages 
and liabilities. 

One such type of signal processing scheme is the discrete 
time or sampled signal processing. This method utilizes a 
Signal which is sampled at regular intervals. The sampled 
values are combined with previous and subsequent samples by 
some algorithm to produce the processed output. Implementa- 
tions of these algorithms are referred to as sampled filters. 
Basic Fourier analysis indicates that periodic sampling in 
time makes the filter response periodic in frequency [1]. 
Busiginies and Dishal report a sampled filter, developed by 
Clark, in 1949 [2]. The repetitious nature of the frequency 
response became useful with the development of radar systems. 

A pulsed radar signal has a discrete spectrum, with 


spectral lines spaced from the carrier frequency at multiples 


eZ 








of the pulse repetition frequency (PRF) [3]. The signal 

to noise ratio of a received signal could be improved by 
passing the signal with the noise through a filter with a 
series of narrow passbands also spaced at the PRF. The 
Signal components would all be passed by the filter while 
only the noise components close to the signal frequencies 
would pass. The remaining noise is attenuated by the filter. 
Such a filter is called a "comb" filter due to the appearance 
of the passbands as teeth of a comb. Figure 1.1 illustrates 


the use of a comb filter. 
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FIGURE 1.1. Signal to Noise Ratio Improvement 
by Use of a Comb Filter 
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While it is theoretically possible to build a comb filter 
out of many individual bandpass filters, maintaining the 
proper harmonic spacing between them is difficult. A 
sampled filter will have a frequency response which is 
periodic at a constant frequency defined as the repetition 
period, f,- A low pass sampled filter, as shown in Figure 
1.2, will be repeated as a bandpass filter centered at 
multiples of fo. hence the sampled low pass filter forms a 


comb filter. 


ee ad ~—omee eee ome 
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FIGURE 1.2. Sampled Low Pass Filter 


For a given continuous low pass filter with a cut-off 
frequency of foe the resulting derived sampled filter will 
have a bandpass bandwidth of 2f.- Since a narrower ‘'tooth' 
will reject more noise, the cut-off frequency of the sampled 
filter should be a small fraction of the repetition period. 
Another application for comb filters is found in radar 


systems which use Doppler shift to detect moving targets. 
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Moving Target Indicator (MTI) radar uses a repetitious 
filter to suppress the ‘clutter’ returns from stationary 
targets while allowing shifted returns from moving targets 
to pass. This canceller filter can be implemented by a 


sampled high pass filter, as illustrated in Figure 1.3. 
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FIGURE 1.3 Sampled High Pass Filter as a Canceller 


The sampled implementation of comb filters may be 
accomplished by a recursive technique where the previous 
inputs and outputs are weighted and fed-back to the input. 


The results of recursive algorithms are treated in great 


LS 





detail by Gold and Rader in Reference 1. The values of 





the weighting coefficients determine the filter cut-off 
frequency. Chapter II of this study looks at the details of 
the relationships between the coefficients and the filter 
response. 

Current literature suggests areas other than radar where 
comb filters are applicable. Arndt suggests the use of a 
comb filters for TV signal enhancement in noisy channels [4]. 
ITlletschko and Mavrodieva utilize comb filters in color TV 
processing of chrominance, luminance, and brightness compo- 
nents {5],{6]. Uses for comb filters in digital controllers 
are suggested by White [7]. Other areas of application will 
no doubt emerge with time. 

Since a comb filter effectively 'samples' the frequency 
spectrum of the input signal with its uniformly spaced 
passbands, information of the signal's spectral content is 
available in each comb tooth. A tunable bandpass filter used 
in conjunction with a comb filter would make it possible to 
isolate the spectral content of a single tooth. Observation 
of the behavior of that signal spectral component in ‘real 
time' would be possible. Webb indicates that cases where 
short time variations in spectral components is of interest 
arise in the aerospace discipline [8]. 

As will be shown in Chapter III, a tunable bandpass 
filter can be implemented by sampled techniques. By taking 
the impulse response of the desired continuous bandpass 


filter response, sampling it at an interval equal to the 
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delay between filter outputs, and setting coefficients for 
weighting the signal samples proportional to the impulse 
response sample values, one may approximate the continuous 
filter response. Such a technique iS non-recursive. Filters 
implemented by the non-recursive method are called transversal 
or finite impulse response (FIR) filters. 

The algorithms used in sampled signal processing tech- 
niques require the operations of delay (storage), multiplica- 
tion (weighting), and addition. The sampled signal values 
may be stored digitally or by analog means. The weighting 
coefficients must also be stored by one of these means. 
Either analog or digital methods may be used to implement 
the arithmetic operations. Until the recent development of 
charge transport devices (CTD), only sampled digital signals 
could be stored. Since digital storage and logic hardware 
has been increasingly available and inexpensive, a great 
deal of study has been done on the digital implementation 
of sampled filters, or simply, digital filters. Now analog 
"shift registers’ are available for implementing sampled 
analog filters. This study wili show that the various 
Z-transform design methods developed for digital filter theory 
may be equally well applied to sampled analog signal proces- 
sors. As a demonstration, a recursive comb filter and a~- 
non-recursive bandpass filter will be implemented and 


cascaded to form an analog Fourier transformer. 


17 








Ii. SAMPLED ANALOG RECURSIVE FILTER 


A. THEORY 

The block diagram of a recursive sampled filter is 
illustrated in Figure 2.1. Since this diagram has two 
Gelays in that particular configuration, it is referred to 


as a second-order recursive filter. 


(<) (+ Bee) 





-b 


+4 \ 


o 
ND 
MN 


FIGURE 2.1 Second-Order Recursive Filter 


18 





Analysis of such a system, performed in Appendix A, indicates 
maeeetrne transfer function, H(z), for this configuration is 


given by, 


-l -2 
. 2 Ves Uae : Ao + a,z + a,z 
i(z) = =—_—_—_—_— = ———_,, 


WeoelZ) -1 -2 
alg! 1 + b,2z + b5z 


where z is defined as Ibe 


and T is the time delay. 
For the special case where a5 = bo = 0, the system becomes 
a first-order filter. Higher order filters may be cbtained 
by paralleling, cascading, or both, first and second-order 
sections. 

The frequency response of a sampled filter can be 


determined from its transfer function, H(z), by using the 


Gemimtion of Zz, 


H(jw) = H(z) | 
z=-e 


ace 

Analysis performed in Appendix B indicates that the delay 
terms introduce sine and cosine variations inthe frequency 
response, H(jw), which yield a periodic nature. The ‘'period' 


of the repetitive response is defined as fos where 


1S) 








Two major design methods are available for determining 
the transfer function coefficients: the standard Z-transform 
and the bilinear Z-transform. Both methods will be studied. 

it otandard 2-Txansform 

The basis for the standard Z-transform comes from 
the impulse response of the desired filter response, h(t). 
The impulse response is sampled at the regular interval, T, 
one time delay, which corresponds to zt in the Z-domain. 
The Z-transform is simply a power series in z+ where the 
coefficients of the series are the sampled impulse response 


values, or 


For many impulse responses of interest this series can be 
expressed in closed form. Tables listing these closed forms 
are called the Z-transform tables. Since the designer often 
knows the Laplace transfer function, H(s), of the desired 
filter, some tables are constructed to allow one to bypass 
the actual impulse response and go directly from H(s) to 
HZ) . 
a. First-Order Low Pass Filter 
For a simple RC low pass filter, the impulse 


response is known, 
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where De aA , the cut-off frequency. This impulse response 


has a Z-transform, 


H(z) = ————_——_ .. 


The 3 dB corner frequency, Wor is found in Appendix to be 


related to the feedback coefficient, bo, by the following 


equation, 
W Ll+b.* - 4b 
on ae -] uF i 
jas 7 a 
s dl 
where We B at Figure 2.2 illustrates this equation. 





FIGURE 2.2 w_vs b 
e Ii 


For a normalized transfer function, ag = 1 and a, = ee 
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b. Second-Order Low Pass Filter 
The impulse response for a second-order LC 


Hew pass filter is, 


h(t) = We sin wot ’ 


-i; 
where Oea= (LC) * . The resulting Z-transform is, 


Wo sin wT zt 
H(z) = SS a ae 
1 - 2cos wt Z. + Z 


For a normalized transfer function, 


ag = Oz; 
a, = os 
a5 = 0, 
by = -2cos wT ; 
bo =1. 


We and by are related by the following equation, which is 


graphed in Figure 2.3, 


ao 








-2.0 -1.6 -1.2 -.8 -0.4 0 0.4 0.8 1.2 1.6 2.0 


by 


BPIGURE: 2.3 Ww, VS by 


c. First-Order High Pass Filter 


The impulse response for an RC high pass filter 


is found in Appendix B to be, 


“wit 

The resulting Z-transform is, 
-wW T 

wee 6 et 
H(z) = = 

= 

Cc =1 
l-e Z 


P48 





By inspection, the coefficients are, 


ag = I, = ey 0 
-w T 
ay = -e S ; 
cn (ed 
ae S 
by = -e : 


Another first-order high pass filter can be 
constructed from the low pass filter by changing the 


polarity of bis 


The polarity change has the effect of shifting the response 
by me so the new corner frequency is = Was where Wo 
and by are related by Figure 2.2. 

d. Second-Order High Pass Filter 


The impulse response for an LC high pass filter 


is given by, 


1, nae 
where Os Te: This yields the following second-order 


Z2-transform, 
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a ie 


1 - (2cos OWT +O. Sin WT) Z + Z 
1 - 2cos wT Z + Z 
thus, 

ao ee ie 

ay = -~2cos wT ~ Wo sin oT ; 

as =il1=, 

by = -2cos WT ; 

Dine = le. 


2. Bilinear Z-Transform 
For the bilinear Z-transform, the Laplace transfer 
function, H{s), is used to determine the Z-transform 
directly by using the following approximation for s as 


described in Reference l, 


-—~ 
N 
I 
} 
nd 





" 
| 
HJ 
N 
4. 
jad 


This formula has the effect of warping the frequency of the 


desired filter response, H(jw), so that, 
= & tan (5) 


iiiees | 
bilinear 


Appendix B shows the precise derivation of this relation. 


Zo 





a. First-Order Low Pass Filter 
The transfer function for a single pole low pass 


filter is, 


Cc Cc -l 
“Eres Tee. — 
e 
HZ) = wT - 2 Qe 
1+ = zt 
wt +2 


A normalized transfer function would have coefficients, 


ay = — 
aT in 

is Lim 2 
b Cc 


1 wT +2 


The relationship between the bilinear Z corner frequency, 


u) and by is given by the following equation, which is 


el 
graphed in Figure 2.4, 


Ipae .). 
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=1.0 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1,0 


dy 


a 8 Os a 8 ee = ee ee 


FIGURE 2.4 w_, vs b 
Cc 1 


Since ag = ay: the bilinear Z low pass filter has a zero for 
the frequency 72W er which means the low pass response goes 
all the way to zero. 

b. Second-Order Low Pass Filter 


The transfer function for an LC low pass filter 


rs, 
Wy 2 
H(s) =—+——— >, 
s+ y 7 
Cc 
= Ty, 
where ia (LC) * . The resulting normalized bilinear 


Z-transform is given by 


_ st 22°72 ae 
Va 

20), 

W 2m + 4 
Cc 


= 
H(s) 5 : 
2 


pe on ) z 9 Ye 


The corner frequency, Wats 1s dependent on by by the following 


relation, graphed in Figure 2.5, 


2a 








eee = i ean}. ene 
W. T 2-b, 
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b, 


FIGURE 2.5 wW_, vS b 
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c. High Pass Filters 
Analysis pexformed in Appendix B shows that 
for the bilinear Z-transform high pass filters, the same 
form and coefficient values with the exception of the 
polarity of a, as in the low pass case are used. The corner 
frequencies are found using the same design curves as for 
the low pass case. 
d. Advantages of Bilinear Z-Transform 
The bilinear Z-transform design method has 
advantages over the standard Z-transform. Given a Laplace 
transfer function, H(s), the Z-transformed transfer function, 
H(z), can be obtained directly. In addition, the bilinear 


Z-transform places a zero to force the stop band attenuation 
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down. The position of this zero may be shifted to transform 
a high pass filter into a iow pass filter, and vice versa, 
by reversing the polarity of the coefficient, ai: The a. 
coefficients of a bilinear Z-transform are always integral 
multiples of one another, reducing deSign complexity. These 
advantages occur at a price. The frequency scaling of the 
bilinear Z-transform response is not linear with the frequency 
response of H(s), but rather folisws a tangent function. As 
a result of this frequency warping, the bilinear Z-transform 
frequency response will only approximate that of H(s) for 
frequencies where, 

Tan — = — : 


U) 
S Ss 


The response cut~off frequency may have to be compensated, 
pre-warped, in order to occur at the desired value. 
3. Sampling Effects 

For either of the design methods described above, 
the frequency response will repeat at multiples of the 
repetition period, fos so that the response shape for 
O<f< af is mirrored from 4f_ < es f and repeated 
from kf. StS (k+1) £. BOI 25 Oys we 9s « FLGUEE TZ. 6 
illustrates the repeated response of an arbitrary filter 


response. 


Ins, 
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Filter 
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af £ 2f 


FIGURE 2.6 Repetitious Frequency Response 


This repetitious response effect is equivalent to considering 
an input at an arbitrary frequency, £, such that 


Pa — |b - Us 5g g8 


1 
< 4f ° 


a S 


That is, the input is '‘heterodyned' with the particular ath 


multiple of the repetition period such that the difference 
frequency, fo: is less than sf jnews1 ter frequency response 
H(f.) is then applied to the input. Therefore, for any input 
Signal frequency, only the filter response from 0 < f. < “sf, 
is used. If the input frequency meets the Nyquist criterion, 


i.e. is less than sf or the multiple of the repetition period 


used to heterodyne the signal, k, is zero. For frequencies 





beyond the Nyquist bandwidth of half the repetition period, 
k is some positive integer. The ambiguity of k is an 
aliasing effect. For k = 0, there is no ambiguity, thus no 


aliasing. 


Bee SAPERIMENT 
f system Configuration 
A first-order recursive filter was implemented using 
a Reticon Corporation Serial Analog Delay (SAD-100) device 
as the delay element. The transfer function coefficients 
were obtained by potentiometer voltage dividers. The 
arithmetic cperations were implemented using operational 


amplifiers with feedback. Figure 2.7 illustrates the system. 
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FIGURE 2.7 First-Order Recursive Filter 
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A variable feedback resistance was installed on the Reticon 

SAD-100 driver card output amplifier so that the gain of 

the delay circuit, -a, could be adjusted. This factor, a, 

is important to the setting cof the coefficient potentiometers. 
Circuit analysis of the coefficient potentiometers, 

performed in Appendix C, indicates that the maximum coeffi- 

client values for ay and by iS or in aaacditron, Since the 

input impedance of the summing circuit is finite, the 

coefficients are determined by a loaded voltage divider 

relationship graphed in Figure 2.8 for the particular loading 


used in the system. 
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FIGURE 2.8 Coefficient Calibration Curve 
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The coefficient calibration curve may be experimentally 
checked by observing the filter response to a unit step 
input. Appendix D derives the expressions for the coeffi- 
cients in terms of the differences between successive filter 
Outputs. Figure 2.9 shows the step response for a typical 
first-order low pass filter. The coefficients are calculated 


by the following relationships, 





This technique is extended to second-order filters by the 


use of the following relationships, 


ag = d, ; 
d.d,d. - dod,” - d.d,d, - 4,74 
_ a) dod, 293 19394 2 ay 
5) 5 ’ 
d,~ - dod, 
ee eee ae ded - ~ d.dad 
ee 24344 2 ds 194 1934, 
2 2 sane | 
d, dod, 
— dyde - dad, 
Dy - 2 : 
a,* - aod, 
2 
ae dy, = dod. 
2 2 : 
d,~ ~ doc, 


Appendix E describes the use of the Tektronix 31 programable 


calculator to solve these equations. 
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For first-order standard Z-transform filters of the 
fee , 


H(z) = —_-+—_, ’ 


1 + bj2 


another method for experimentally determining by 1S available. 
As discussed in Appendix B, by taking the ratio between the 
response minimum and maximum, U, bi is determined by, 


|= 
1 Lltu f 





|b 


which 1S graphed in Figure 2.9. 
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2. Experimental Peequency. Response 


Several types of filters were implemented to verify 
the design theory and observe the characteristics of tne 
repetitious response. For an implementation using a delay 
time, T, theory indicates that the filter response should 


repeat exactly at multiples of fo: where, 


Since the SAD-100 is a 98-bit register, the delay time is 


98 times the clock period, Toe 


Ee oO Tae 
e 


The repetition period would then be the clock frequency, 

oe divided by 98. Since a new input is sampled every clock 
cycle, the rate of sampling is fa: The sampling is done at 
the input gate of the SAD where a transistor is turned on 

by one of the two clock phases for half of its period. 

Since the two phase clock operates at af os the ‘on' time, 
i.e. the sampling pulse width, is To: The finite pulse 


width causes an envelope over the sampled signal spectrum. 


The envelope follows a =="— shape. This envelope has its 
fest null at f = = = — »- Since the clock frequency is 
Cc 


N times the repetition period, where N is the number of 


delay bits, N frequency response cycles will occur before 
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the first sampling envelope null. Because of aliasing from 
the sampling rate, foe only those responses for frequencies 
less than or equal to fs iver eiewrirst oN, are useful; 
In this study, only the first twelve cycles. are observed, 


therefore the minimum value of the envelope was, 


sin n (2% 
98 


Soa! 


which is sufficiently close to unity to be neglected. 


Standard Z-transform filters were implemented with, 


—OMeo — =dese-0./ Low pass, 
by = 025,07 O.7 © ALoOW pass. 


Bilinear Z-transform filters were implemented with, 


a_=l, a, =1, b, = -0.7, 0.0, 0.3, 0.7 low pass, 


=-l, b, = -0.3, -0.7 high pass. 
Two clock frequencies were used, 

f = 400 kHz, 1 MHz . 
Reoures 2.10a through 2.19a illustrate the first response 


period of the various filters implemented. The solid curve 


represents the calculated response while the points are 
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the observed response. Figures 2.10b through 2.19c illus- 
trate several resvonse cycles of the filters. 

Some Significant deviations from the theoretical 
responses were noticed, especially for the larger values of 
b,- The responses are not mirror images of one another. 

The reason for this 'skewness' is not clear at present. One 
possible explanation may center on a frequency varying 
coefficient, such as aj The values of ay and a, must stay 
equal for a zero to appear at multiples of f.- If one of 
the values changes differently with frequency than the 
other, say through a frequency characteristic of the delay 
device, the zero would not be properly formed, perhaps 
shifted up or down frequency from kf... The other major 
deviation is the change of the maximum amplitudes of the 
responses with frequency rather than ree uniform as 
theory predicts. 

For filters with small b,, the trend is for decreasing 
amplitude with increasing frequency, particulary beyond 40 kHz. 
This effect is probably due to the gain roll-off of the 
operational amplifiers in the summers. Part of this decrease 
could be due to the sampling envelope, though it is always 
@eeaeer than 0.97. 

For bilinear Z filters with by = +0.7, the maxima 
of the frequency responses increase with frequency, as 
shown in Figures 2.11 and 2.14. This trend may also be due 


to deterioration of the zero from ay and ay varying more 





with increased frequency. The effects of frequency varying 
coefficients need to be pursued. 

Measurements of corner frequencies, wnich are 
summarized in Table I, indicated good agreement with the 


design theory. 
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FIGURE 2.10a Theoretical and Observed Frequency Response 
ef Standard Z Low Pass Filter; by = -0.3 
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FIGURE 2.10c 
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FIGURE 2.lla Theoretical and Observed Frequency Response 
of Standard Z Low Pass Filter; by = -0.7 
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FIGURE 2.l1lb Theoretical and Observed Frequency Response 
of Standard Z Low Pass Filter; by (0) 37 
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of Standard Z High Pass Filter; by = 0.3 
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FIGURE 2.12b Theoretical and Observed Frequency Response 
of Standard Z High Pass Pilter; Dy = 0.3 
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of Standard Z High Pass Filter; by = 0.3 
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of Standard Z High Pass Filter; by = 0.7 
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FIGURE 2.14a Theoretical and Observed Frequency Response 
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FIGURE 2.14b Theoretical and Observed Frequency Response 
of Bilinear Z Low Pass Filter; by = -0.7 
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FIGURE 2.15a Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter; b, = 0.0 
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cf Bilinear Low Pass Filter; by = 0.0 
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FIGURE 2.16a Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter; b, = Ora 
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FIGURE 2.16b Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter; by = 0.3 
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FIGURE 2.17a Theoretical and Observed Frequency Response 
Of Bilinear Low Pass Filter; b, = 0.7 
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FIGURE 2.17b Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter; ws) 
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FIGURE 2.17c 
of Bilinear Low Pass Filter; by = 0.7 
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FIGURE 2.18a Theoretical and Observed Frequency Response 


of Bilinear Z High Pass Filter; bj = -0.3 
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FIGURE 2.18b Theoretical and Observed Frequency Response 
of Bilinear Z High Pass Filter; Dy =O 
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FIGURE 2.18c Theoretical and Observed Frequency Response 


of Bilinear Z High Pass Filter; b, = -0.3 
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FIGURE 2.19a Theoretical and Observed Frequency Response 


of Bilinear Z High Pass Filter; by = -0.7 
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FIGURE 2.19c 7 
of Bilinear Z High Pass Filter; b, = -0.7 
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Filter 


Bilinear Z 


Low Pass 
by = 0.0 
0.3 
0.7 
High Pass 
by = -0.3 
-0.7 
Standard Z 
Low Pass 
by = -0.3 
-0.7 
High Pass 
by = 0.3 


TABLE I 


Theoretical 


elo y, 


0.056 


Qa22 5 


O70 57 


0.279 


0.443 


OeZo 
02358 


0.438 


Oe boS 


OR0GL 


DAS, 


0.072 


0.254 


0.429 


O85 IE Sy 


02065 


ORAL, 


OPA OKS 6 


G.252 


0.428 





IIIT. SAMPLED ANALOG NON-RECURSIVE FILTER 


Awe) LHEORY 
Since the frequency response of a sampled filter is 


periodic, it may be described by a Fourier series [ll], 


co 


H(jw) = 2 D. ej UTk ; 
k=-0 
where 
1/T Sie 
D. = — f mG) Go a 
-7/T 


This expression for D. is mathematically equivalent to an 
inverse Laplace transform of the transfer function, H(s), 


which is h(t), evaluted at t = kT, 


Di = h(kT) , 


where h(t) is the impulse response of the filter. By 


using the definition of Zz, 


Hezjye =. F h(kT) me ’ 
k=-e 


where h(t) is assumed to be an even function in time, 


therefore it is not a true impulse response. 
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For a delay device with 2N taps available, H(z) can be 
approximated by the series, 
N 


Hitz) = A BUKT) =& 
=-N 


k 


Figure 3.1 illustrates the implementation of this series. 
The approximation can be made closer to the designed 
response by truncating the series at larger values of N. 
At least 2N storage elements are required. 

For a rectangular bandpass filter response, such as 


that shown in Figure 3.2, repeated with a period of 


ea 
We 7 f 


lated in Appendix G by the following expression, 


the coefficients of the Fourier series are calcu- 


7 a. 0 
D,. = Diy are Sin(2kTBT) cos(2k ra : 
Since truncating the series is like multiplying the 
Mieeonee Series by a window which is umaty for -NT < t < NT 
and zero elsewhere, the result of the truncation is to 


convolve the frequency response with the Fourier transform 


of the window function. A rectangular window transforms to 
Sin x 

x 
a result of the convolution. By choosing a different window 


function which adds sidelobes to the response as 


function, such as a Hamming window, these sidelobes can be 
Minimized at the expense of greater passband width [8]. To 
use a Hamming window, each coefficient, Dy is muitiplied 


by a weighting factor, Whe where, 
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FIGURE 3.1 Implementation of a Non-Recursive Filter 


Ue 





H(¥) 


2B 
“WoW -W -W 4+W -W TE v4 
8 Oo 8 3 0 16] Oo oO W<-W S W+wW 
s oO s 


FIGURE 3.2 Rectangular Bandpass Filter Response 


W. = 0.54 + 0.46 cos(n =) 


kK 
Since the input to the non-recursive filter is sampled 

with a finite width pulse of T,, the clock period, the 

spectrum of the input is enclosed within an envelope, (aes 


sin(7f£T,) 


E(£) = Te ° 
C 





This envelope has its first null at £ = fae tne clock 
frequency. Though this envelope is imposed by sampling 
rather than the filter algorithm, it can be included within 


the filter response function as a multiplicative factor. 


(ie 





Eee GAPERIMENT 
ie COnLIGuUration 
Reticon Corporation markets a 12-term Tapped Analog 
Device (TAD-12) which provides 12 output taps with the 


following time delay configuration. 





FIGURE 3.3 TAD-12 Delay Configuration 


Since the first tap is only ave delayed from the input and 
the remainder of the taps are spaced with 2T ue the undelayed 
input can not be used for the A, term, therefore, only ll 
terms of a given tapped delay device can be utilized. 

The filter is implemented by weighting the tap 
outputs proportional to the Fourier series coefficients. 
Analysis performed in Appendix G and the computer program 
described in Appendix H provide design values for the 
tapping resistances, R, Figure 3.4 describes the filter 
set-up. The lower ends of Ry are connected either to the 
inverting or non-inverting input of the summer depending 


upon the relative coefficient polarity. Since D. = Diy s 
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FICURE 3.4 Non-Recursive Filter Set-Up 
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R, = R.,- %$In order to not load down the taps, a minimum 
value of 10,000 ohms is set for the tapping resistances. 
Appendix H indicates the following expression for Rye 


R, = (lox + x) (88%) -xr, 
Dy. | 
where r is the tap output impedance. For the TAD-12, r is 
about 5K. Since the resistances depend on the relative 
values of the series coefficients rather than the absolute 
values, gain scaling is accomplished by adjustment of the 
feedback resistance, Re - For consistent scaling for 


al 


Opposite polarity coefficients, R must be equal to R 


f as 

Due to the TAD-12 EPeces the tap outputs are at 
a +5 volt bias level, so the lower end of = returns to 
+5 volts rather than ground. The output of the summer, 
therefore, rides on a5 volt bias level. 

2. Observations 

Two sets of coefficients were designed for a Q of 

15. The first set are for fo = “af and the other for 


= kf . Table II lists the coefficients and corresponding 


: S 


0 
resistances. These sets of coefficients were measured at 
clock frequencies of 50 and 100 kHz. 

Figures 3.5 through 3.12 illustrate the experimental 
frequency response compared to the theoretical calculations. 


For these data, the measured center frequency, foe passband 


width, BW, and Q are listed in Table III. 
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fy = \ 
|k | Dik 
0 Ones 333 
iL 0.0 
2 =).027/70 
3 0.0 
4 Une Os 55 
5 Gr 0 
fo PES 
oO Z 0.25 
0.5 
100 kHz Oe 5 
0.5 


TABLE II 


Ss o ~ ts 
SK | [| SIX | 
10 K C2 03533 10 K 
open -0.03084 eZ 
Py a.0 0.02419 i lie ee 
open -0.01556 Ze) cl 
85.0 0200773 Don / 
open -0.00314 LS4°5.3 
TABLE III 

18 BW 

O 

5.25 KHZ Gee 3. KHz. 

LZ S2063 
BA 6.402 
220 6.154 
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FIGURE 3.5 Theoretical and Observed Frequency Response; 
i = 4f 0 ; fo = 50 kHz 
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FIGURE 3.7 Theoretical and Observed Frequency Response; 


Ea = 4f. ’ fo = LOM als 
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FIGURE 3.9 Theoretical and Observed Frequency Response; 


fy = ae ‘ fo = 50 kHz 
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FIGURE 3.10 Theoretical and Observed Frequency Response; 


f, = af ; £ = 50 kHz, Magnitude in dB 


8 3 





1.0 


0.8 


= 
ON 


Magnitude (normalized) 
Se 
f» 


0.2 





\ 
pea eee a a eas 
20 40 40 50 





Frequency, kHz 


FIGURE 3.11 Theoretical and Observed Frequency Response; 


iG = aie ; Eo = 100 kHz 
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FIGURE 3.12 Theoretical and Observed Frequency Response; 


Lo _ af. : fe = 100 kHz, Magnitude in dB 
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It is noteworthy that Qdid not change significantly when 
the clock frequency was doubled and that the bandwidth 
changed by a small amount when the center frequency was 
doubled by coefficient changes. 

The measured Q was Significantly below the designed 
value of 15. Some insight into the reason for this discrep- 
ancy may be found by looking at the impulse response of a 
rectangular bandpass filter. Without rejecting the negative 
time portion of the impulse response, since the non-recursive 
filter uses a ‘two-dised' impulse response, one has from the 
inverse Fourier transform of a rectangular window centerea 


at plus and minus Wor 


W) sSin(t — 


= cos (wt) [ 


2 
Q w 
S 


This is a continuous function for a non-repeating frequency 
response. For a frequency response periodic with period, 


W the impulse response is discrete with the above function 


oo 
acting aS an envelope. Figure 3.13 illustrates the response. 
In actuality, the sampled nature of the impulse response 
causes the repetitious frequency response, but the duality 


of the Fourier transform ignores the true cause and effect 


relationship. The spacing of samples is, 


De 


E [NN 
IS 
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FIGURE 3.13 Two-Sided Impulse Response 
W 
For — = a , there will be a 1 samples per cosine period. 
S 


There are Q divided by pi cosine periods under the main lobe 
of the envelope. Therefore, the number of samples within 


the main lobe is, 


If one uses the criterion that the main lobe must be 
represented completely in order to approximate the impulse 
response, i.e. M samples are required, then some estimate 


of the number of terms required to give a selected Q at a 
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center frequency of aw, 1s available from the above 


equation. Conversely, if the number of terms is fixed, the 
realizable Q may be estimated given the ratio of W, to wa. 
In this study, M= 11 and a= 0.25 and 0.5. Table IV 


lists the expected and measured Q's. 


TABLE IV 
wo 4 5 
W expected measured 
OFZ5 8.64 oS 
OS ees 4.1 


The observed Q's did occur in about the expected ratio of 
two to one, but were significantly less than the expected 
values. Analysis in Appendix J indicates that this lowering 
of Q is mostly due to the Hamming weighting used. The 
expected values for this weighting are 2.17 for the center 
meecmency tO repetition period ratio of 0.25 and 4.35 for 
the 0.5 ratio. These match more closely to the observed 
values. The remaining difference is probably due to the 
need for more impulse response terms than those in the main 
lobe. 

A significant amount of clock noise, up to 40 milli- 
volts, was observed at the output. A low pass filter would 


be desirable at the output to decrease this noise content. 
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Perhaps this low pass filter could be built in conjunction 
with a level shifter to remove the output bias. It was 
further observed that using lead lengths of a few inches 

to connect the coefficient potentiometers to the TAD added 
Significantly to the clock pick-up, which severely limits 
the small signal handling capability of the device. Input 
Signals beyond 8 volts peak to peak began causing distortion 
in the tap outputs. Therefore, without using a low pass 
filter to suppress the clock noise, the dynamic range of 


this device is less than 54 QB. 
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IV. SAMPLED ANALOG FOURIER TRANSFORMS 


fee THEORY 

As mentioned in Chapter I, a possible application of 
sampled analog signal processing is to use the combination 
of a recursive comb filter and a non-recursive bandpass 
filter as a sampled analog Fourier transformer. The teeth 
of the comb filter breaks a signal into its Fourier compo- 
nents, while the bandpass filter selects one of the 


components. Figure 4.1 illustrates such a svstem. 








Comb Filter Bandpass Filter 
H(f) 
Out 
ai 
S(f) 
’ \ 
/ > as A 
f : » f Lf \ fo \." f 
Input Spectrum Sampled Spectrum Selected Tooth 


FIGURE 4.1 Sampled Analog Fourier Transformer 
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B. EXPERIMENT 

For the purpose of demonstration, a bilinear Z-transform 
low pass filter with a narrow tooth width and a non-recursive 
bandpass filter are implemented and cascaded. For this 


example, the comb filter coefficients are chosen as, 


ay = 1, 
ay = 1 ’ 
by = -0.7 e 


The design of a generalized analog Fourier transformer 
can be made based on the number of frequency components 
desired, i.e. the number of comb teeth selected. The number 
of desired components, K, determines the frequency of the 


highest order tooth to be selected, 


Max S 


Since ee = =, where T = LT. for L equal to the number 
of delays in the recursive filter delay device and ae is 


the recursive filter clock period , 


Max 


Since f. is the sampling rate, the Nyquist criteria states 


that, 


sdk 





Max 


So 


For a bandpass filter designed for a center frequency to 


repetition period ratio, 


u) 
2 — a 
a 
aS 
the next highest passband is located at (1 - a) w. co 2S uenee 


the first K teeth should be distinct, the first repeated 


passband should occur at least at a frequency, 


Ke ye 
while, 
Ay 
AW. ae 7 ° 


These restrictions on a can be summarized, 


b+ 





If the bandpass is tuned by changing its clock frequency, 


iecomO wile ronan concetant. SO its bandwidth will increase 


o2 





as it is tuned to higher ordered teeth. At the maximum 


value, i.e. K f , the bandwidth will be, 





The term, 8, is a design choice that should be some fraction 


less than unity. The required Q is, therefore, 


K 
Re ai ‘ 


From the equation relating Q, a, and M, the number of taps, 


a guess at the required device length is made, 


K (K + 2) 


Q 
a = aaa : 


& TI 


Thus, for a possible design to find 100 frequency components, 


and a maximum bandwidth of one tooth repetition period, 


z.e. B=Hal1l, 
M > 200(200%2) 5947 
This means that the Fourier series is from the ~1623°? term 


Perthe 1623°° term. 


The non-recursive filter coefficients are selected for 
a center passband frequency of af Figure 4.2 shows the 
resulting transfer function for the cascaded filters plotted 


against the comb response. The comb filter clock was run at 


a3 





1.0— 


FIGURE 4.2 
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at 400 kHz, which placed comb teeth at approximately 
4.1 kHz increments. The bandpass filter clock was adjusted 
so that the passband selected the third tooth. The frequency 
of the third tooth maximum was measured to be 12.452 kHz. 
The bandpass filter clock was measured to be 49.555 kHz. 
Several aspects of the cascaded response are noteworthy. 
The shape of the passband is determined by the comb filter 
tooth shape, which should remain constant from one tcoth to 
the next, whereas the non-recursive bandpass filter bandwidth, 
as shown in Chapter III, changes as it is tuned with its 
clock. This use of the comb filter to maintain a constant 
spectrum ‘sample' window allows one to keep higher order 
harmonics separated. The other important aspect of the 
system which is noteworthy is the isolation achieved between 
the selected tooth and the unselected teeth. The following 


table lists the first several teeth and the isolation of each. 


TABLE V 
TOOER Isolaticn 


mcon OB 

se aa rane 
0.0 

2 Ore 7 

sot 

-44.0 

= 8) If oll 

4s Oo 

soo We 

-44.0 


NH 


Owoann u& Ww 


bw 


The low isolation figure for the ninth tooth resulted from 


the repeated response at 5 fe : 


a> 





V. COMCLUSTIONS 


Gne Z-transtorm Cesign methods developed: for digitel 
filter theory are useful for sampled analog recursive filter 
design. High and low pass cut-off frequencies were predict- 
ably determined by the filter feedback coefficient values. 


These values were adjusted by means of a generalized cali- 


bration curve or by measurement of the unit step ressonse 
of the filter. Some deviations from theoretical frequency 


responses were observed. These deviations are apparently 
caused by instrumentation limitations. 
Non-recursive filters were implemented by truncation of 


a Fourier series expression for the desired repetitive 


Oy 
(D 
af 
y 
Os 
UO 
~ 
i 


frequency response. The center frequency of t! 


filter response passband was tuned at nearly constant Q by 


tty 


Varying the clock frequency. The pass band was also shifted 
at nearly constant bandwidth by chances mace to the coeffii- 
erents of the filter. The sin X over X envelope caused by 

the finite sampling pulse-width multiplies the desired 
frequency response. This result decreases the sicnals passing 
through the repeated passbands. 

Several particular frequency components of an input signal 
spectrum were examined in real time by cascading a recursive 
eeomo. filter and a non-rectirsive batdpass f1i¥%ter to form 
sampled analog Fourier transformer. An estimate of the 


~ 


a sampled 


Cr 
tty 


required delay device lensths for implementaticn o 


Fourier transformer by the above method was found for a genera 


number of desired components. 


le 


= 


Bi 





APPENDIX A 


Derivation of Filter Transfer Function 


For the system diagram illustrated below, 





FIGURE A-l 


the transfer function, H(z), is defined as, 


rE ee 


a Seen | 


Ha: 
From the system diagram it follows that, 


=2 


X(z) = B, (z) + (-by) (X(z) 2-7) + (-by) (X(z) 2), 
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SO, 
peeecZ)) = X (Zee X(z)z > + b X(z)z2- 
an 1 2 
The output of the system is, 


E out '2) = AX (z) + A,X(z)z 7 + ayX(z)z * x 


therefore, 


7 Eee Cz) - X(z) (ag + az > + Anz ) 
CE) NS Se ae | ee 
in X(z) (1 + biz 


ay = 
Ay oT a,z a= aoa 


1 + biz 2 
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APPENDIX 8B 


Recursive Filter Design 


The transfer function of a given second-order sampled 


recurSive filter is given by the following expression, 


Peete ae 
0 1 "2 
A) ees 2 a 
Je ap D,2 + b.z 
By uSing the definition of z, 
Z= aps = eIJut , where T is the time delay, 


the frequency response, H(jw), can be determined: 


-jwT = 7 ie 
ay + a,e + aod 


1 + b,e 


zZ 


ag + a, (cos wl -j} Sin wT) + a, (cos 2wT - j Sin 2wT) 


1 +b, (cos wl - 3 sin wf) + bj (cos 2wT - 3 sin 2ut) © 


Since the only variation of H(jw) with frequency is due to 


the sine and cosine functions, and since they are periodic 


in w for a fixed T, H(jw) is periodic in frequency. [In fact, 
the period of the response is (k) at < w < (k+l) aa ' 


we, 





mee K = 2ee-2,-1,0,1,2,-+-- « Since the sine and cosine 

terms repeat themselves exactly, the entire frequency response 
of the filter can be described by looking at a single period, 
say 0 < w <a. The frequency, a , 1s called the 
repetition period, Wer which is related to the time delay, T. 
As illustrated in Figure 2.6, the frequency response from 

0 <w< %W, is mirrored for 40, <4 < Ww, . This periodic 
Nature of the filter frequency response, which arises from 

the use of delays, can be used to implement comb filters as 
described above. 

There are several methods for designing the coefficients 
of the transfer function, H(z). Two methods are frequently 
used if the Laplace transfer function, H(s), or the impulse 
response, h(t), of the desired continuous response is known. 
One is the standard Z-transform procedure. The other is the 
bilinear Z-transform method. 

For the standard Z-transform approach, h(t) is determined 


for the desired response and then sampled at intervals, T, 


the time delay, 


See 


D sampled 7 nigh) OG Cie) : 


lor 8 


0 


Its Laplace tranform is, 


skT 


2 cone i mee A({kT) ¢ ; 
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From the definition of z= ac ; 


H(z) = £ h(kt) 2%. 


Z-transform tables are available for those impulse functions 
of interest where a closed form of the series is known. 

The second technique, the bilinear Z-transform, involves 
the use of an approximation for s in terms of z. The 
transfer function, H(s), 1s converted directly to H(z) by 


the substitution of, 


N 
! 
}—! 





49) 
i 
le | Eg 
N 
be 
a 


The use of this approximation for s leads to a non-linear 











relationship between the frequencies of s = jw and 
re | 
oe ee a 
aoe. 2 2a 1 
Jw Toe 
SO 
d -}; cei —j44w'T 
J wt zZ(z*-2z *) _e -e 
EL Sie a —ahy'T ’ 
@ 2*(z* +2 2) atau TT, Jaw t 
therefore, 
Se = tan ‘w'T , 
Or, 
2 1 Se tes ae 
oS tan (4w' T) and ow! = — tan “(4oT) . 
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eal 


Wes 
wy’ = 7 tan (OL) ae 


This non-linearity is defined as frequency warping. 


A. STANDARD Z~-TRANSFORM DESIGN PROCEDURE 
l. First-Order 
a. Low Pass Filter 
A first order RC low pass filter has the impulse 
response, 


= eoke 1 
ett) € ; where w= Ra 


The standard Z-transform is, 


Il 
= 
mM 
D 
N 


H (z) 


Ag = WO 2 
a, =O, 

‘ ~w it 
by = -e : 


IL9)% 





Since the value of ay 1s a constant factor which sets the 
filter gain, normalizing the response is equivalent to making 
a normalized transfer function with ay = 1. Since the 
coefficient by 1s the only term which acts on the frequency 
G@ependent z, the cut-off frequency, Woe is determined solely 


by b,- From the definition of cut-off frequency, 


& 


ilk : 
2 1H(3) |oax 


, 2 
[H(jw,) | 


1 
1 + b, (cos wT = ap eae WT) 


1 


cos wT + b é 


it 2 1 


He 


The maximum value of H(jw) occurs when 2b, cos WT 1s 
most negative, i.e. that factor is equal to ~2|b,| . 


Therefore, 


i = + 
Zn 2 
1 + 2b, cos wT + b, een |e 2 
Or; 

1 + 2b, coswT+b.% =2- Align let 2 é 

di E iE L ul 
os ee ae 

a . : 

COS WO. Ib 


1 


naOks: 





2 i) 
ee = We = cos ~ ( 2B, 

i, 29415.) BS ee 
cee cos” + ( | 1 = 
W. 27 2b, 

TABLE VI 
s 
S10) 7 Us 0.5000 
-0.2 0.3524 
=O 0.2207 
-0.4 O2nS 73 
= 55 Of EES 0 
=O 0 020631 
=e 0.0574 
10) 5 ts: Oe s57 
=-0.9 0.0168 
=e) 0.0000 


For the first order low pass filter, by a0 ene re heme 
maximum values will occur whenever 2b, cos wi = 2b, PE 858 


when cos wl = 1, or 


D: 
Ww = k (=) >, Jeo SE a ee 


Kw, r 
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Likewise, the minimum value of the response occurs when 


wT = -1 t or ’ 


e 
it 


(2k+1) (5) 


(k+%2) w 


The ratio of the maximum response to minimum response yields 


an interesting result, 


db 
> (Hse) |[%,, 1 + 2|b,| + by’ 
a MH (3u) |2 a I 
J max 2 


1 - 2|b, | + by 


(lo! [ay Ns 
cain TP ee 
(1 + |b, |) 
thas |b, | 


Solving this equation for b yields, 


- 
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TABLE VII 


u |b, | 


0.0 1.0000 
0.01 0.9802 
O5025 029512 
0.05 0.9048 
0.1 0.8182 
OeZ 0.6667 
0.3 OP sy Sitei5) 
0.4 0.4286 
OB 023325 
0.6 0.2500 
O/ 0.1765 
0.8 Se eden 
0.9 0.0526 
1.0 0.0000 


This result can be of use in experimental determination 
of the effective value of by: 

b. High Pass Filter 

The low pass filter just described above can 

be transformed into a high pass filter by shifting the 
response by 73W so that the maximum values occur at 
(k + 45) w. and the minimum values at Ku. Changing the 
polarity of by would accomplish this result since the 
maximum response occurs when 2b, cos wT = -2|[b,| , which 


would make cos wf = -l1 and w= OS ee) Lor by > 0 
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The minimum response would occur when cos wl = 1, or 


= Kw. The cut-off frequency, W ae is found by the same 


procedure as for the low pass case: 


Z 
max 


JH(ju) 1* = 518 (34) | 


(cere en Ee 
“jw ot 


1 tbe 


IL 


2 
cos WT + b 


io) 2b 1 


HN 


From the previous case, Since the expressions for the 


cut-off magnitudes are identical, 


ee ee fs - 4b, + a, 
W 27 2b, 
TABLE VIII 
are 
b, “s 
0.1716 0.0000 
0.2 0.1476 
0.3 0.2793 
0.4 0.3427 
0.5 0.3850 
0.6 0.4169 
Oh 0.4426 
0.8 0.4643 
0.9 0.4832 
1.0 0.5000 


OY 





HOorTtnisS feommeot the first order high pass filter, the other 

coefficients are the same as for the low pass case, i.e. 

ag = 1 a, = 0 for a normalized transfer function. 
Another type of high pass filter may be imple- 


mented by the standard Z-transform technique. For a first 


order RC high pass filter the impulse response is, 





h(t) =f tH(s)] =L tS 
Cc 


cles 
e 7 C = 
=—emattt) Woe ; where w, = 


The Z-transform of this impulse response is, 


ee ewe es 
H(z) = xX (6(nT) = wee ) z 
n=0 
-y T 4 
i- We ae C 2 1 
= are T ~ = 
l1-e Z 1 


The transfer coefficients are, 


eo) 
ec: 

eo Be ’ 
-w _T 

b, = -e : A 


The corner frequency is a function of both ay and bj. 
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2. Second-Order 
a. Low Pass Filter 
For a second-order LC low pass filter, the 


impulse response is 


2 
u) 
h(t) =o ie (=I = W, sin wt , where ie (LC) 2 


Cc 


The resulting Z-transform is, 


co 
= : me 
HZ) =") 2 Ww, sin w nt Z 
n=0 
w sin wT zt 
aoe C 
iSe2ecos (ie zits, 
For a normalized transfer function, ag = ag = OF 4 by =e 
and, 
a, = a 
by = -2 cos wT ;: 


From the expression above, the relationship between by and 


the corner frequency, Was can be stated, 


ete 


ag 
= gee OS 1) 


IONE, 








TABLE IX 


2534) 0.0000 
oy Se! 0.0718 
=At6 0.1024 
alee 0.1266 
0.1476 
0.1667 
late 0.1845 
OR PONIES: 
Occ 
OG 5 BAL 
Oh ASIC, 
052659 
e220 
OR2905 
O21 55 
OR 3556 
0.3524 
0.4734 
0.3976 
0.4282 
025000 


' ot 
a eae 
e oe 
© No 


ee 
NORM HHH eH OOo CO oe © © 
@ 


COON RPNOAAVHN HPN ON ff OV 


b. High Pass Filter 


The impulse response for a second order LC high 


pass filter is, 


2 
h(t) = £7 [+S] = 6(t) - ow sin wot . 
S une 


This transforms into, 


2) re) (S§(nT) - We Sim nw T) zg 2 
n=0 
-i -2 
7 1 - (2 cos WT + w sin WT) Z + Z 
iL -2 


1 - 2 cos wT 7A 
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where, 


ag = i 

a, = -2 cos WT + Wo sin wT 
a, = Ly 

by =Eaer COs wT ; 

bo es 


The cut-off frequency for this filter will be dependent on 


both ay and by. 


Be BILINEAR Z-TRANSFORM DESIGN PROCEDURE 
l. First-Order 
a. Low Pass Filter 
The transfer function of a first-order RC low 


pass filter is, 





N 
j 
a 





il 
Hi] bt 


By substitution of the bilinear approximation, s 


N 
a 


lil 





U) 


Hl 


RZ) 


kr} 


Lo 
Za 





2 
= + 
T Wo 


fe 


w T(z +1) 
2Z - 2 + wz +1) 


= 
Wot +e WoT Z 


_ = 1h 
2t+ Wo + Cw, ZZ 


wT . WT vl 
2+, 2+, 
WT - 2 


ee -1 
wT +2 








1 + Z 


ed: 
Normalizing the transfer function by ae yields, 
fe. 





= 
+ 
H(z) = —+ 4 — ’ 
n -1 
ics = a, 
Ayo eee: 
where, 
ag = Li 
= ie Ls 
wy T-2 
b c 


1 wt + 2 


The bilinear Z-transform warps the frequency response, so 


UNewactUal Cut—-OoLf frequency, Wor is , 
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From the expression for Dy in terms of Was 


ee 





SO, 
Oa ee Ib 
Wg: = tan GD) 
CL, 
mc = 7i tan | 
W TI Del : 
S 1 
TABLE X 
Wat 
25} =e 
Ss 
-1.0 0.0000 
= 058 O-0167 
~Oec 020352 
=a 0.0556 
=Ora6 02.0780 
S055 0.1024 
-0.4 O212389 
ood UP) Oro 2 
AOA ONS 2 
—I0) 5 Jt O23 3 
0.0 U6 PASLONG, 
Os OZ oily 
Qe 2 Ore eas: 
O73 0.3428 
0.4 Oras aL 
O85 Gime SiS) JAS 
OG 0.4220 
Om? 0.4444 
0.8 0.4648 
Oa 0743833 
0 0.5000 


JE 16S 





b. High Pass Filter 
For a first-order RC high pass filter the 


transfer function is, 





Ss 
H(s) = ’ 
STW. 


which transforms to, 





Ea 
<P 
im eo 7 Me 
7 z-l1 
tyes 


For the normalized transfer function, 


ag = 1 a 

a, = -~l] , 
— Wit 2 

1 wit +2 : 
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It is noteworthy that the only change from a high pass filter 
to a low pass filter, other than the gain normalization 
factor, is the polarity of aj: Since the expression for by 
in terms of the design cut-off frequency, Was 1s the same 
as for the low pass fiiter, the same relationship between 
by and the actual cut-off frequency, Wor as found for the 
low pass case exists. 
2. second-Order 
a. Low Pass Filter 
The transfer function for a second-order LC low 
pass filter is, 
w 4 7 
H(s) =——“-— , where fe = (Cues Saar 


2 
So +0 
Cc 


mae bilinear Z-transform for this function is, 


2 


Wo 


H(z) ————— 
22-1,2 Z 
le Sa! FW, 





we T (2°+z41) 
2 22 


4zZ° -~- 8z + 4 + We (oo 42 241) 








where, for a normalized transfer function, 


ay = il. & 
= i 
a, = Ae 
2 Bet - 8 
b = 
ak 22 44 
Cc 
bo =]. 


By using the expression for by to solve for 2 Inte eine Os 


bi, one finds, 


ma p ha see 
@ T 2-b, 


Because of the frequency warping, 


DES eS 


Mees a oe 
Mets oe ao a ieee! ae 
Or 
Wa Z iL eae sat A 
Q) T 2-b 
S lL 


It is interesting that this design table is the same as 
Table IX, the standard Z-transform second-order low pass 


filter. 
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TABLE XI 


b “c' 
ih W 
-2.0 0.0000 
-1.8 0.0718 
~1.6 0.1024 
-1.4 0.1266 
=e. 2 0.1476 
~1.0 0.1667 
-0.8 0.1845 
-0.6 0.2015 
-~0.4 0.2180 
-0.2 0.2341 
0.0 0.2500 
0.2 022659 
0.4 0.2820 
0.6 0.2985 
0.8 0.355 
1.0 GO233Ss5 
es 2 0.3524 
1.4 0.3734 
1.6 6.3976 
1.8 0.4282 
2.0 0.5000 


b. High Pass Filter 
A second-order LC high pass filter has the 


following transfer function, 


2 
H(s) ==>, 
Ss +W 
Cc 


which is transformed to, 
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Fic) = Tzttl 
22-1,2 2 
le S47! as are 





Ae eo aay 


42° - 824440 °T* (2742241) 


4 E:) ie ie 
en? 44 Seas ar ad 
= c Cc Cc 
2% 2 ae 
1 + a) Z, + 2 
aT +4 


For a normalized transfer function, 


ay = l=, 
a, = 2.77 
as = d Nie 
2W 2m? _ 
pee ; 
= WW) 2n2 44 
Cc 
bo =l. 


It 1S noteworthy that this transfer function has the same 
coefficients as the low pass filter with the exception of 
the sign of ay: The corner frequency is dependent on Dy 
in the same manner as for the low pass case. It is now 


apparent that the polarity of ay is the means for converting 
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high pass bilinear Z-transforms to low pass and vice versa. 
In addition, the Agr ays and a, values are always integral 


multiples of one another. 


Cc; SUMMARY 
The following table summarizes the design procedures 
for the standard Z-transform and the bilinear Z-transform 


for first order filters. 
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Design Method 


Standard Z 


Bilinear Z 


Filter Type 


TABLE XII 











Coefficients 
ay ay Dy 
“yw T 
WO 0 -aaic 
wT wot 
l-w -e e 
Cc 
-—w T —w) T 
Cc Cc 
1 -w -e -e 
Cc 
wt 
Ww) 0 e 
Cc 
wt wt - wt -2 
w T+ 2 w T+ 2 WwW T+2 
Cc Cc Cc 
wt OT w T= 2 
wT + 2 wT +2 wi t+ 2 
2 Pe hee 
+ 
wot + 2 wT + 2 wT 2) 
, 9 wT -2 
+2 
wt + 2 ut + 2 wt 
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Period 














APPENDIX C 


Derivation of Calibration Curve 
For Recursive Filter Coefficients 


The coefficients of the recursive filter transfer 
function are produced by potentiometers. These potentiometers 
are loaded by the input impedance of the summing circuit and 
are driven by the output of the SAD-100 driving board as 


illustrated below. 


Vin = R 


(Ears R 
—— 
R, Vout 
19 


FIGURE C-l Coefficient Circuit 


It is assumed that the summer op amp is operated as a unity 
gain, inverting device. Therefore, the voltage on the tap 


of the pot is -v - The voltage on the top of the poten- 


out 


tiometer is ON tte The coefficient, b, is defined as the 


matl@ OL eae to Vane From the voltage divider equation, 


E21 





a's = SO. [RFR TT ? 
out in R, + Ro Fe 





therefore, 
aR R 
7 Bee R,+R 
a -V 7 R_R 
in Ry ns a 
R, +R 
RR 


The limit of the coefficient, b, as R, the loading resistance, 


gets arbitrarily large is, 


R,R Ro 
limb = lim 6 Sea oO 
Roe ie ae eee: Ro 


This is the expected value for a non-loaded voltage divider. 
Since R is finite in this study, its effect must be 
considered. For the actual hardware used, a, = 10K, 

R= 30.9K, anda = 1. The following calibration curve was 


calculated, 


122 





rH oO 
A 


mo ON HD WI & W NH 


_ 
© 


0.000 
0.097 
0.190 
0.281 
O2n i 
0.463 
Osa) 
0.655 
0.761 
0.875 
1.000 


TABLE XIiil 


123 


epeerved 


0.000 
02075 
0.180 
0.286 
0.369 
0.454 
0.546 
Or 030 
0.734 
0.858 
1.000 





APPENDIX D 


DETERMINATION OF TRANSFER FUNCTION 
COEPPIeCIENTS ROM UNIT © STEP RESPONSE 


A particular filter configuration can be characterized 


by a transfer function, H(z), where 


a, ta Z 
0 1 
EeCZ) = iA ToT, ae 


1 + by Zz + b 


sT 


for z defined as e and T a time delay. Since a transfer 


function relates the output of a system to its input, 


E (z) 
H(z) = -ey , 
Stopes 
Or, 
Fut 2) = Ef?) H (z) 


-1 _ 
ay a5 Baie + ae 
e) e 


—] a 
1+ biz + bp 2 


= Bint?) ° 


Multiplying the above equation by the denominator of H(z) 


yields, 


~] -—2 _ 
E out (2? * by E out 2) Z + bo E out 62) Z = 


ie 


-] 6 
ao EB. (2) + a) Es, (2) z, + a5 Ey 2) Zs; 
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SO, 


Taking the inverse Z-transform yields the following time 


expression, 
Cour (7) = Apes y(n) + a,e,,(n-l) + aoe. (n-2) = bye oe (BoD) 
— boeoye (nS) 
where, 


e (i) = 


out (iT) . 


“out 
Assuming zero inputs and outputs for n < 0 and a unit step 
mimpoutc, 1.€. e.,() = 1 for all positive n, the output of 
the filter will be the following sequence for the first 


five values of n, 
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Table xIV 


e577) Sous 

0 0 0 

1 1 Ao 

Z il ay + a, - bia 

3 ir ay + ay ta, - b, (apta,—-b,a))- boay 

4 1 a ta,tas—b, (apta,ta,—b, (apta,—-b, a) 
~b5a,))-bs (a ta -b, a) 

5 au ajta,ta,—b, (apta,ta,—b, (aptaj ta, 


~b, (ata ,—b,a,)-b,a.-b, (a ta,—b, a5) ) 


~b, (agta +a 


1 97b, (ata -b_a,)-b 


17189) “P24! 


Defining the difference between successive steps in the 
output wavefore as d = Gout $2) = Cour tinh. the following 


relations are found, 


d. = a.ta.ta 


gt ta2-b, (ayta,—b, ap) —bya,-ayta,—b, a, 


ay ~ b, (a,-b, a) - ba, : 
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Similarly, 


dy = -b, (ab, (a, —b) a9) -b5 a9) - b, (a,-b,a)) 


The relationship for higher differences is, 


i, 2 op Cl a= By Chi Oe Soe 


It is clear that the tern, ayta,ta,, 1s subtracted out in 
each difference calculation after k = 4, thus the higher 
differences are a function of only bys bos and previous 


differences. Solving this relation for by yields, 


= Crake 1G 


by = eee P for k > a 
ads 
SO, 
a oe eo ee iy 
i d, 


From the higher difference relation, 


den —OAdS 


Solving Lor bo and substituting the above relation for bys 


one finds, 


2) 





Gem = dd 


sehr ene 
3 24 
Substituting this expression into the relation for by 
yields, 
ae So eer 
i 2 
oor 9 204 


The a coefficients are now easily determined, 


- dod, — 4,44, 
Ee 
oe eee owed 
3 224 
Beil Bi GeGle RUBLE ee! 
ons 293 193%, ey 
2 t 
dd, - dy 9, 
ao = d., + bid, + bod, 
ie cecmcnets ae aid ~ - a.d.a 
Sexe 293%, a Ss SIS eee 
5 
fy och ey 


All coefficients are now determined solely as functions of 
the differences between successive steps in the output. 


Higher order differences could be used to determine by 
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and bo, but using dq, and d 


5 


minimizes the differences required 


eondetermine the five coefficients. 


In the case of a first-order transfer function, 


ae. 


a, = bo = 0, a simplified set of equations would result, 
ay = dy ; 
a, = d 1 bj 4d, 17 18) dy 
0 = d., + b,d, + O-d, ; 
SO, 
d 
3 
b. =l-—_— r) 
il d. 
a, = dq, + bid, 
2 
wee os 
2 
Furthermore, since dd = “bd, 4 “bd, _5 and bo = 0, 
ee 
ce 


This is a useful result for experimental determination of 


the transfer function coefficients through measurement of 


the step response since by may be determined for several 


values of k and then averaged to reduce error. 
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APPENDIX E 


PROGRAMS TO CALCULATE THE RECURSIVE FILTER 
CORE PWetENTS EROUSTHE UNIT STEP RESPONSE 


The formulae resulting from the derivation in Appendix 
D lend themselves to a programmable algorithm for solution. 
The programmable Tektronix 31 calculator is used to imple- 
ment the algorithm in the laboratory. The first five step 
differences from a unit step input are entered in order. 
The calculator then determines and prints the five second- 
order coefficients of the transfer function, Agr Ayr Ags 


b eana-b.. 


td 2 
The basic instructions for using the program are listed 
below: 

1. Turn on calculator and press “clear"™. 

2. Insert tape cartridge containing filter 
coefficient programs. 

3. Press "reset" and transfer the desired program by 
pressing "from tape" and the number of the track 
the program is stored on (0 for second-order, 
Z2e5Or first-order) « 

4. When "busy" light extinguishes, enter d, on 
keyboard and press "cont". (Note: the calculator 


Seoul sorrnt Out CG. -. ) 


1° 
5. Repeat step 4 for dq through d. The calculator 


Ve Veeriniohiwexecutson Of the program after “cont” 


B54) 








is pressed for d (di, for first-order 


=) 
program). 

After the coefficients are calculated and printed, the 

program returns to its beginning and is ready to read a 

new d,- 

Two limitations of the second order program are note- 
worthy. First, the determination of all coefficients 
except A is dependent on bo which is, in turn, dependent 
on the higher order differences. For the coefficients of 
interest in this study the differences became steadily 
smaller, therefore, the higher ordered ones contain larger 
percentage error. This error then propagates through the 
other coefficients. Careful measurement of the step 
response can minimize this problem. 


The other limitation of the second-order program concerns 


the factor a, - dod,- For a first-order filter, 


3 4 
b=-r-=~F=r- =>, 
1 d. d, 
Or, 
see 
Ga = 4495 
thus, 
7 = 
d. = d., qc, = Q. 


ich 








This factor, however, appears in the denominator of the 
expressions for all the coefficients except €@g: Even ag 
experimental error kept the denominator from equaling zero, 
the answer would be useless. For this reason the first- 
Order program waS written. As noted above the first-order 
program requires four differences to be entered. The 
Soecrtficient, bys is calculated by the following weighted 


average, 


The average weighting is a crude compensation for the 


Meeater percentage error in the d, measurement. 
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VDE PDO 1S 


COMPUTER GENERATION OF THEORETICAL 
RECURSI Ve VETETER FREQUENCY RESPONSE 


For a given recursive filter transfer function, H(z), 


where, 


2 a a 
l + by Z + bo Z 


the frequency response, H(jw), may be calculated by using 


the definition, 
z=e =e = cos (wl) + j sin (wT) , 


where T is the time delay. By Substitution, 


a, + a, (cos wl - 3 Sin wT) + a, (cos Zul =} Sin 2h} 


0 
a b, (cos Qi e je sin. wl act b. (cos 2QL toedeslnecul) 


H (jw) 


tid=cos Ml + acos ZuT = j(a,sin wl + a,sin 2wT) 


2 2 
cos 2wT - j (b,sin dy a b.sin 2uwT) 


a0) 1 


1+ b.cos wt + b 


I Zz 


For a given frequency, the magnitude and phase of the filter 


response can be determined, 


2 
= [ReH(jw)}]* + [ImH(jw)]° , 





| H(3a) 


pee tink (9) W) 


£H(jw) = Tan [Rel (ja)! : 


iS 3 





A computer program was written to calculate, given the 
transfer function coefficients, the response magnitude and 
phase over a specified number of steps of a specified 
interval of frequency. The results were then plotted on 
the Calcomp plotter. Experimental magnitude versus fre- 
quency measurements were read in, frequency normalized, and 
super-imposed over the theoretical curves. 


A sample data deck follows: 


1. Line Printer Plot Characters 


ee 





UZ 8e a 8S 


2. Transfer Function Coefficients 
J r.0 -1.0 > 0. -0.7 ot ENO) 
/ 


| 
2... WU... 2021... 3031... 40 41. . 50 
| 


All numbers F10.5, in order: Agr ay, 


3. Experimental Data Head Card 


v4 9 0 4.125 


/ 





2 Soa yes 13 


nS JAM Os: 
# Of points Response Repetition Period (kHz) 
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4. Experimental Data Points 


. Pe aa = wens a oe a a ee ee ee ee 


Ue 3557 0.001 


| 


co citi GIL fee 7410. 


: beOne BAO eS 
Frequency (kHz) Magnitude (normalized) 


The program will print on the line printer the coeffi- 
cients, the data points, a theoretical magnitude versus 
frequency plot, a theoretical phase versus frequence plot, 
and a message stating the status of the Calcomp plot, if 
requested. A labeling statement is used to title the 
Calcomp plot with the coefficients of the transfer function. 
The title declaration cards need to be changed as different 
coefficients are used for the title to coincide with the 


response determining coefficients. 
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APPENDIX G 


Non-Recursive Filter Design 


The design of: a non-recursive filter is based upon the 
periodic nature of a sampled filter frequency response. 


The response may be described by a Fourier series, 


ey 
Jk (==) w 
H(jw) = 2 Dy e ' 
k=-co 
where 
uy Aik 
oe — esate ane 
k 27 
qe es 


The particular filter response of interest in this study 1s 


a bandpass filter with rectangular passpands, 


H(f) 


- i er nae 


FIGURE G-l Rectangular Passband Filter 
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which can be described with the parameters: response 
repetition period, fi passband center frequency, Eo? 
and pass bandwidth, BP. The filter quality factor, Q, 


is defined as, 


Q= pp! 


and the fraction of the response repetition period which 


corresponds to the center frequency is defined as, 

FFS = : 
The filter description is complete by specifying fo. Q, and 
FFS since the center frequency and bandwidth can be calcu- 
lated from tne above three parameters. The repetition 
period is, 


fe =i FZ OG ate radians/second. 


For a rectangular bandpass filter, 


iL Fae | ee 


0 | 
H(jw) = 
0 elsewhere - 


therefore, 


bf 





2m (£5+%BP) 


a: . go Sku? 
De 27 J : ‘ os 
2m (£)-*BP) 
> ~jk(2nf£5T) 
= io sin(kt Be sip SC 
T 
= <- sin(kn BP T) (cos 2nf T = 3) label 2nf£,T) : 


Since only real terms may be implemented, the series coeffi- 


cients are, 


2 
k i sin(k7n BP T) cos (21f)T) 


D = 
The series expansion has infinite terms; however, the non- 


recursive filter has a finite number of taps available to 


approximate the series, 


I 
tJ 
\e) 


H (jw) 


™1 
Se) 
* 
a 
O 
9) 
x 
E 
Z 
UJ. 
9) 
t’- 
5 
x 
E 
= 
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N 
H(z) = z D, Zz : 
k=-N x 


Since the Z-transform must contain only non-positive powers 


of z, i.e. no anticipated signal values only delayed values, 


N 
lal e4) 1 2 N = » D. Zeon z 


k=-N 
This relation states that the transfer function will contain 
a constant delay of NT from input to output in addition to 
the response phase. A possible implementation of this functicn 
is illustrated in Figure 3.1. From the expression derived 
for the bandpass filter, 
Ce ora lL Sin (ky BP.) 


7 Mr eecT csr itieer “eo Co™ Or’ 


it can be seen that this is an even function of k, so, 


In addition, 


lim Dy, =e bee tu l)  (1)) =.2:BP-T , 
I atz 


this can be expressed in terms of the design parameters, Q 
and FFS, 


D = 2°BP’T =2--—~-: 


ie9 





Window weighting functions are used to modify the 
coefficients in order to improve the filter characteristics. 
A Hamming window is frequently used where the weighted 


coefficient, Ayr HES, 
A, = D. (0.54 + 0.46 cos (* T) ) x 


The expression multiplying D. is called the window weighting 


Bunction. 
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APPENDIX H 
Computer Calculation of Coefficients and 


Tapping Resistances for the Non-Recursive 
Filter 


From the analysis performed in Appendix G for a bandpass 


feLeer, 
= = FFS sin(kxX) _ 
Dy = Diy = 2 a 0b cos(kY) , N< k<N, 
where, 
BP a 
X = 7 I BP = pass bandwidth 
Ss f. = response repetition period, 
EG 
Y = 27 Fo fo = passband center frequency, 
Ss 
= Eo 
Q = Ep ’ 
fe 
= Se 
FFS = a : 
S 


For a Hamming weighting window, 


A, =A 


k 
re = D,. (0.54 + 0.46 cos (— 1) ) : 


-k 


For the TAD-12 analog delay device, the coefficients are 
implemented by input resistances to the operational amplifier 


summing circuit, as shown in Figure H-l. 
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Vv 
out 





FIGURE H-l Coefficient Implementation 


In order to minimize the loading effect on the tap, the 
tapping resistors, Rye are greater than or equal to 10 K ohms. 
mae contribution of a single tap, Vir to the output, ree 


ise Che product of an and the coefficient, A,, or, 





k Ve ames R, 


This expression indicates that R. is inversely proportional 
to A|, therefore the largest coefficient will correspond to 


the minimum allowable resistance, 10K, so, 
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Solving this equation for R, and substituting that result 


into the Ay expression yields, 


A (r + 10K) 


_ max 
=== po 
k 
Or, 
R, = ——= (r + 10K) Soe 
JA, | 


For the TAD-12, the output impedance for a tap, r, is about 
SK. 


A computer program was written to calculate the weighted 


and unweighted coefficients and their respective tapping 


resistances for a given T = aes Pp Opto, anata palr Of N. 


These parameters are inserted directly into the program as 


follows, 











YW 0.25 
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4. Truncation points, NM1 and NM2 
Fa 2 
1 eo ¢ 6 7 s 9s s 
/ NM2 = 11> 
1 es 2. ¢ 6 7 s s se 
Note: The program calculates results for, 


NM1 


H(z) = £ A, 7 K-NM1 
k=-NM1 
and, 
NM2 
H(z) = » A, 2X MN z 
k=-NM2 


5. Maximum passband magnitude, CA 


Note: CA was set to one in the derivation in Appendix 
G, but the program is more general in this 


respect. 
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APPENDIX I 


Computer Calculation of Non-Recursive Filter 
Frequency Response 


Since the transfer function for a non-recursive filter 


is, 
N 
H(z) = Ay 2 , where z=e ; 


it 1s possible to calculate the magnitude and phase of 


H(jw), where, 


N 


H(jw) = © AY J SUH 
k=-N 
N 
= =) A, (cos (kwT) + 7 eSin (ur) ) 
k=-N 
N 
= Ay + oy [(A, +A_,) cos kwT + J (A, -A_,) Sin kwT J 


The real part of H(jw) is clearly the cosine terms, 


Re H(jw) = 


in S&S 


Ap + fas (A, +A_,) cos(kwT) , 


while the imaginary part of H(jw) is the sine term portion, 


N 
Im H(jw) = 2 (A, -A_}) Sin(kwT) . 
k=1 
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The magnitude and phase of H(j ) are readily calculated, 


lH(jw)| = [Re H(jw)]* + [Im (jw) ]? 


1 


jH(ju) = tan”? (eR HOY) 


RS Eee 


A program was written to calculate the magnitude and 
2 


phase of H(jw) from O = jw = = , plot the result on the 
line printer and Calcomp lotter for a given set of A,- In 
addition, experimental data points may be read and plotted 
over the theoretical curves. 

The above analysis assumes perfect impulse sampling, 
however, actual devices will have some sampling envelope, 
usually a Sin function, which will multiply the magnitude 
of the response. Provision is made in the program to 
introduce this envelope. 

The following sample deck will illustrate the use of the 
program. | 

imeeetnc Printer Plot. Characters 
a ee Y 


l32345 


§ 
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2. NMAX, the number of coefficients 


NMAX = 2N + 1 


eae 
(I3) 


3. Coefficients 


(F10.5) 


4. Data Header Card 
NPTS, the number of data points (multiple of 30) 
FCK, the clock frequency 


XX , the value of sampling envelope at 
magnitude maximum frequency 


[oe 
/ 6 0 and) 0.9845 


iit 14 ay 4.23 


NPTS BCK XX 
(I3) CELLO 5) EEO 25) 
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DE bata! Olles 
FX, observed frequency 
HX, normalized magnitude 
aa 6.2735 0.033 


es kl ee 20 


FX HX 
| (F10.5) (F10.5) 


Note: If dummy cards are required to make NPTS 


a multiple of 30, use FX = 999., HX = 1. 


When the Calcomp plotter is used the "call draw" cards 
are included and the title declaration cards are modified 
to title the plot with the proper number of terms and FFS 
value. The program will always print the number of coeffi- 
clients, the coefficients, and the data points so that proper 


data entry can be checked. 
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APPENDIX J 


Effects of Hamming Weighting 


The use of a weighting factor, Whe in the time domain 
is equivalent to convolving the frequency response of the 
desired filter with the Fourier transform of the weighting 


function. The Hamming weighting function, 
= =e 
Wy (t) = (0.54 + 0.46 cos (sa ke) ; 


which has a Fourier transform, W, (3) ; 


G 


W 
« pen ° * S ® ieee S 
W, (34) = 0.54 é(jw) + 0.23(6 (jw + 55) + 6(jwW ay) s 


where, 


U) _ 27 
Ss ‘i a 
This result can be illustrated as three impulses as in 


Figure J-l. 


0.54, H(w) 
0,23 0.23 
-W 0 W % 
eae. 8 
2N oN 


FIGURE J-l Fourier Transform of Hamming Function 
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When a bandpass filter 1S convolved with the impulses, its 
bandwidth is increased by an amount proportional to the 
spacing between the impulses. This added bandwidth will be 
a fraction, y, of the Hamming function bandwidth, 


Ys “Bs Ys 
added BW = Y(5y 7 ) = OL oe i 


2N 
The value of gamma is estimated by the following procedure. 
If the values of the impulses were equal, then the -3 dB 
points of the convolved response would be shifted, roughly, 


W) 
by spacing of the impulses, = , aS illustrated below. 


H(w) 


—_— = aod — . 
. 


-W -W O- w+Ww W 
Cc crs 


FIGURE J-2 Equal Impulses 


The added bandwidth would be equal to the spacing between 


W 
S 
the impulses, Wo Chea 
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For the case where the outboard impulses are zero amplitude, 
the response is just convolved with a single impulse and is 
exactly reproduced, so zero bandwidth is added. In this 
case, gamma is zero. A relationship is suggested by this 
bracketing of gamma, 

amplitude of outboard impulses 


tome amplitude of inboard impulse : 


enus £Or a Hamming function, 


y= == = 0.426 


Since the Hamming weighting adds some to the bandwidth, B, 


the Q is expected to be less, 


ea a. 
BW + B ° 


By dividing the expression for Q' by that for Q yields 


iD Sal 





BW 


Q' = 2 ays 3 


BW 
YW 
elie Sipe set 


aN 


an + yO ’ where a = 


ae 


al os 
ro) 


9) 


It is noteworthy that the expression for Q' indicates that 
the limit of Q' as N gets larger is Q. 

If one uses the main lobe criterion for expressing the 
impulse response by specifying the responses within the first 


nulls, one finds that, 


Q = (2N+1) amt 2 


Substituting this result into the expression for Q' yields, 


(QN* +N)q7 


| — 
he (2m7y +1)N + yr 


For the case of N = 5 and y = 0.43, 


Q' =a (8.69) : 


Semtouear— O0e25.and 0.5, 0" will be 2.17 and 4.35, respectively. 
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